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SUMMARY 


An  investigation  is  made  of  nonlinear  oscillations  in  which  the  damping 
and  static  moments  are  represented  by  arbitrary  polynomial  functions  of  the 
dependent  variable .  When  the  nonlinear  damping  is  small  but  the  static  nonlin¬ 
earities  arbitrarily  large,  an  approximate  solution  is  established  which  leads 
to  expressions  for  the  damping  decrement  involving  elliptic  integrals  and 
gamma  functions  in  special  cases.  An  "effective  linear  damping"  is  defined  and 
a  generalized  formula  for  this  parameter  is  obtained  that  is  valid  for  a  wide 
range  of  nonlinearities  in  both  the  damping  and  static  moments.  This  formula 
is  useful,  for  instance,  in  deducing  the  dynamic -stability  parameters  of 
missiles  observed  in  nearly  planar  motion  in  free  flight. 


INTRODUCTION 


Many  physical  systems  are  describable  only  in  terms  of  nonlinear  ordinary 
differential  equations.  One  example  is  the  large -amplitude  pitching  motion  of 
hypersonic  aerodynamic  configurations  for  which  the  frequency  can  be  strongly 
dependent  on  amplitude.1  In  the  analysis  of  ballistic-range  data,  it  is 
desired  to  deduce  the  aerodynamic  properties,  linear  and  nonlinear,  from  obser¬ 
vations  of  a  given  set  of  oscillations.  Toward  this  end,  a  knowledge  of  the 
effects  of  the  nonlinearities  involved  in  the  differential  equations  governing 
the  motion  is  very  important. 

In  this  study,  we  will  be  concerned  chiefly  with  the  damping  of  oscilla¬ 
tions  of  nonlinear  systems.  Although  we  will  be  interested  in  solutions  to  the 
nonlinear  equations  for  their  own  sake,  the  eventual  goal  will  be  to  derive 
formulas  allowing  the  extraction  of  both  linear  and  nonlinear  damping  param¬ 
eters  from  a  set  of  data.  From  the  standpoint  of  aerodynamics,  this  goal  is 
equivalent  to  determining  from  data  the  dynamic- stability  parameters  of  a  mis¬ 
sile  oscillating  in  a  plane.  A  theoretical  study  indicating  the  form  of  the 


^Although  the  present  analysis  is  applicable  to  a  wide  variety  of  oscil¬ 
latory  phenomena,  the  text  is  cast  in  terms  of  pitching  oscillations  of  bodies 
having  nonlinear  damping  and  static  moments  in  flight  in  a  medium  of  constant 
density. 


nonlinearities  to  be  expected  in  the  dynamic -stability  parameters  can  be  found 
in  reference  1.  Aspects  of  obtaining  free -flight  damping  in  a  conventional 
wind  tunnel  are  discussed  in  reference  2. 

Since  we  will  be  interested  in  problems  in  which  some  of  the  nonlineari¬ 
ties  may  be  large ,  conventional  methods  of  analysis  in  which  all  nonlinearities 
are  considered  small  (see,  e.g.,  refs.  3  and  4)  are  not  always  suitable.  Two 
methods  that  treat  large  nonlinearities  in  the  static  moment  have  been  devised 
by  Murphy:  a  perturbation  technique  (refs.  5  and  6)  and  a  quasi-linear  tech¬ 
nique  (ref.  7)  that  yields  results  essentially  the  same  as  those  of  Coakley, 
la it one ,  and  Maas  (ref.  8).  Although  these  methods  are  quite  general  in  that 
they  are  applicable  for  combined  pitching  and  yawing  motions,  only  cubic  static 
moments  are  discussed,  and  explicit  damping  formulas  for  only  linear  damping  in 
planar  and  near-circular  oscillations  are  derived  and  plotted.  Redd,  Olsen, 
and  Barton  (ref.  9)  derived  formulas  for  the  effective  damping  of  nonlinear 
oscillations,  but  considered  the  static  moment  to  be  linear. 

Many  aerodynamic  systems  are  not  describable  by  cubic  static  moments,  and 
quintic  or  higher  order  polynomials  are  often  essential  for  proper  data  reduc¬ 
tion.  There  is  reason  to  expect  that  if  the  static  moment  requires  description 
by  higher  order  polynomials,  then  the  danping  moment  will  also.  It  is  desir¬ 
able,  therefore,  to  understand  the  influence  of  arbitrary  nonlinear  damping  and 
static  moments  on  oscillatory  motions.  The  intent  of  the  present  paper  is  to 
contribute  to  this  understanding. 

The  analysis  proceeds  along  two  distinct  lines  after  appropriate  transfor¬ 
mations  are  used  to  render  the  original  second-order  differential  equation 
into  an  equivalent  integral  equation  on  the  one  hand  and  an  equivalent  first- 
order  differential  equation  on  the  other.  By  means  of  the  integral  equation, 
for  which  a  special  exact  solution  exists,  an  approximate  solution  is  found 
that  is  valid  for  large  danping  nonlinearities  and  small  effective  static  non- 
linearities.  For  small  damping  but  arbitrarily  large  nonlinearities  in  the 
static  moment,  an  approximate  solution  is  then  found  by  utilizing  the  first- 
order  differential  equation.  For  this  latter  line  of  analysis,  the  damping 
decrement  can  be  conveniently  represented  by  an  "effective  linear  damping," 
which  can  be  evaluated  in  terms  of  quadratures  leading  to  elliptic  integrals 
and  gamma  functions  in  special  cases.  Based  on  limiting  values  of  the  quadra¬ 
ture  expressions,  a  simplified  approximate  formula  for  the  effective  linear 
damping  is  devised  that  is  extremely  accurate  for  a  wide  range  of  nonlinear 
parameters  and,  further,  is  amenable  to  data-reduction  techniques. 


NOMENCIATURE 


A,B  constants  of  integration 

a  amplitude  variable  defined  by  equation  (28b) 

b  defined  by  equation  (8lb) 

D0  defined  by  equation  (57c) 
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E 

f 

H 

Hoe 

Hs 

Hs 

h 

hs 

I 

K 

k 

M 

Ms 

Ms 

ms 

P 

Q 

R 

s 

t 

to 

u,z 


complete  elliptic  integral  of  second  kind 
arbitrary  function  occurring  in  equation  (27) 


s>o 


effective  linear  damping  defined  by  equation  (66) 

damping-moment  coefficients 

defined  by  equation  (ll) 

H0/v 

i  H. 


rT  yQ  (eq.  (54)) 


2  +  s  HQ  “'a 
defined  by  equation  (8lc) 
complete  elliptic  integral  of  first  kind 
modulus  of  elliptic  integrals 
L„  s+i 

V 

s>o 


Hn 

M  -  — 


static -moment  coefficients 


Ms  -  Hs 

MsymS/v2  or  MsyaS/v2  (eqs.  (42),  (5*0) 
defined  by  equation  (l8) 
defined  by  equation  (30a) 
defined  by  equation  (3*0 
summation  index 

independent  variable  in  equation  (l) 
initial  value  of  t 

functions  defined  by  equations  (3),  (7) 9  and  (8) 
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y  y2/ya2 

YQ  value  of  Y  for  zero  damping 

Ym  mth  order  expansion  term  in  Y  (eq.  (56)) 

y  dependent  variable  and  amplitude  of  nonlinear  oscillation 


v 

T 


a  particular  maximum  amplitude 

nth  maximum  amplitude 

(n+l)th  maximum  amplitude 

initial  amplitude  of  y 

zeroth  order  approximation  (eq.  (50)) 

defined  by  equation  (79b) 

tan“1(H0/2v)  or  sin'^H 0/%J%) 

m2(2m0  +  m2) 

\iAo\l/2  if  Mo  ±  0,  |Mnyn|l/2  if  Mq  =  0 
frequency  variable  defined  by  equation  (28c) 


co  angular  frequency 


BASIC  ANALYSIS 


The  Basic  Equation 

In  this  study  we  wish  to  consider  oscillations  that  are  influenced  by  non- 
linearities  in  both  the  damping  and  static  forces  or  moments.  We  shall  assume 
that  the  nonlinearities  are  functions  of  the  dependent  variable  only. 


Consider  a  nonlinear  oscillation  to  be  governed  by  the  following  differ¬ 
ential  equation  for  y  =  y(t) : 
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The  middle  term  represents  the  damping  moment,  the  last  term,  the  static 
restoring  moment.  We  shall  refer  to  the  term  represented  hy  s  =  0  in  both  the 
damping  and  static  moments  as  the  linear  moment.  Similarly,  we  shall  refer  to 
polynomials  terminated  by  s  =  1,  s  =  2,  and  s  =  4  as  quadratic,  cubic,  and 
quintic  moments,  respectively.  When  Hs  and  Ms  are  nonzero  only  for  even 
values  of  s,  the  motion  is  called  symmetric;  when  Hs  and  Ms  exist  for  odd 
values  of  s,  the  motion  is  called  unsymmetric.  Symmetric  motions  are  of  the 
greater  interest,  but  unsymmetric  motions,  such  as  oscillations  about  a  trim 
angle  in  aerodynamics,  do  occur. 

Equation  (l)  includes  as  special  cases  several  well-known  equations.  If 
H0  i  0,  Mq  i  but  Hs  =  Ms  =  0  for  s  ^  0,  equation  (l)  is  linear  and  has  a 
well-known  solution.  If  H0  =  -H2  =  -e  and  Hs  =  0,  s  ^  0,  2;  Mq^O  and 
Ms  =  0,  s  ^  0,  then  equation  (l)  reduces  to  the  classical  Van  der  Pol  equation 

^  +  e(y2  -  1)  §  +  -  0  (2) 

at2  a 

If  all  the  damping  terms  vanish,  leaving  only  a  nonlinear  restoring  moment, 
then  the  motion  is  periodic.  If,  in  addition,  the  restoring  moment,  the  last 
term  in  equation  (l),  is  a  quintic  or  less  in  y,  then  exact  solutions  may  he 
found  in  terms  of  elliptic  functions  for  symmetric  oscillations. 


Equivalent  Pair  of  Equations 

It  is  convenient  to  seek  a  transformation  that  will  express  equation  (l) 
in  such  a  form  that  the  first  derivative  of  the  dependent  variable  does  not 
appear . 

Let  us  introduce  two  functions  u(t)  and  z(t)  so  that 

y(t )  =  u(t)z(t)  (3) 


Since  we  have  introduced  two  functions  to  replace  one,  we  have  an  arbitrary 
condition  at  our  disposal.  We  will  select  u(t)  so  that  the  differential 
equation  for  z(t)  does  not  explicitly  involve  its  first  derivative. 


If  the  transformation  given  by  equation  (3)  is  substituted  into 
equation  (l) ,  then 


dfz  +/  2  du  + 


dt' 


u  dt 


HsuSzS 


dz 

dt 


s>o 


i  afs  +  VYm  +  i  uV 

u  dt2  \  s  u 

s>o 


z  =  0 


(4) 


It  does  not  suffice  to  set  the  coefficient  of  dz/dt  in  equation  (4)  equal  to 
zero  in  order  to  specify  the  needed  equation  for  u.  If  this  were  done,  a 
term  involving  dz/dt  would  appear  in  the  coefficient  of  z  owing  to  the 
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presence  of  d2u/dt2.  The  form  of  the  coefficient  of  dz/dt,  however,  suggests 
the  following  relation  for  u: 


IdU  1  \  „  S  S  /  r-\ 

u  dt  “  “  2  as^su  z  (5) 

s>° 

Here  the  as  are  constants  to  he  determined  so  that  the  dz/dt  terms  vanish 
in  equation  (4) . 

Taking  the  derivative  of  equation  (5)  to  determine  d2u/dt2  and  substitut¬ 
ing  into  equation  (4),  we  find 


z  =  0 


(6) 

We  now  eliminate  the  derivative  dz/dt  from  equation  (6)  by  requiring  that  its 
coefficient  vanish.  This  can  be  done  term  by  term  if 

_  2 

as  2  +  s 

Hence, 


dfz 

dt2 


/'(-as+l- 2  aB) 


Hsuszs 


s>o 


dz 

dt 


s>o 


Ms  + 


a  c*  *"* 


A£s  du 
+  1)  u  dt 


s  s 
u  z 


dfz 

dt2 


+ 


2 


fs _  1  du\ 

+  s  u  dt ) 


uszs 


z 


=  0 


(7) 


(s> 

s>o 

Equations  (7)  and  (8)  are  the  ones  we  have  sought.  These  equations 
together  with  equation  (3)  replace  the  single  equation  (l)  for  y  by  two 
simultaneous  equations  in  the  new  variables  u  and  z.  It  may  appear  that  the 
solution  of  the  two  simultaneous  equations  for  u  and  z  would  be  more  compli¬ 
cated  than  the  single  equation  for  y.  We  will  find,  however,  that  the  equa¬ 
tions  become  uncoupled  for  special  values  of  Mg  and  Hg,  and  exact  solutions 
may  then  be  found.  Moreover,  transformations  of  these  equations  will  yield  a 
single  integral  equation  for  y  and  a  single  first-order  equation  for  y  that 
are  appropriate  for  establishing  approximate  solutions. 

To  see  more  clearly  the  nature  of  equations  (7)  and  (8) ,  let  us  write 
them  in  a  slightly  different  form.  Substituting  equation  (8)  into  equation  (7) 
gives 
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z  =  0 


(9) 


d2z 

dt2 


'( 1 


H  \2 

s  uszs 


2  +  s 


s>o  \s>° 

Equation  (9)  may  be  simplified  by  introducing  new  parameters,  Hg,  defined  by 


I*-v  ’( l 


H  \  2 

s  s  s' 
u  z 


2  +  s 


(10) 


s>o 


s>o 


The  square  of  the  series  in  equation  (lO)  may  he  arranged  so  that  Hs  may  he 
identified  as 


H 


s 


HjnHs- 


m 


s  /  (2  +  m)(2  +  s  -  m) 

m=o 


(11) 


Equation  (9)  becomes 


d2z 

dt2 


^(Mg  -  Hs)uszs 


s>o 


z  =  0 


(12) 


With  equation  (7)  written  in  this  form,  special  exact  solutions  suggest 
themselves . 


Special  Exact  Solutions 
It  can  be  seen  that  when 

equation  (l2)  reduces  to 


where 


Ms  =  Hs,  s^O 

(13) 

to 

0 +  (Mo  -  t)z  ■  0 

(14) 

solution 

z  =  cos  T 

(15) 

T  -  V(t  -  tQ)  ,  V2  =  Mq  -  ° 

7 

and  tQ  is  an  arbitrary  constant.  The  other  arbitrary  constant  which  would 
normally  appear  has  been  suppressed  since  it  will  appear  in  the  solution  of  the 
equation  for  u(t) ,  equation  (8) . 


Consider  now  a  two-term  damping  moment  given  by 


Here  we  have  kept  the  linear  term  HQ  and  one  nonlinear  term  of  arbitrary 
power  n,  where  n  is  an  integer.  With  this  damping  moment,  equation  (8) 
takes  the  form 


1  du  _  Hn  ,,n„n 

u  dt  '  "  2  "  2  +  n  U 


(17) 


A  solution  to  equation  (17)  may  be  found  by  introducing  the  transformation 


U2  = 


-Ht 
P  e  ° 


(18) 


so  that 


i-  B 


2  dP 


dt  2  +  n 


(19) 


which  has  the  solution 


H„t 


P=U2e°  =  I  A  + 


nH. 


n 


2  +  n 


zn  e 


dt 


(20) 


where  A  is  a  constant  of  integration 


From  equations  (10)  and  (13)  we  find  that  for  this  special  case  the  only 
nonzero  values  of  Ms  and  Hs  are  given  by 


M„  =  Hn  = 


H0Hn 


n  n  2  +  n 


(21) 


H. 


Man  =  H2n 


n 


(2  +  n) 


Finally,  by  combining  the  solutions  for  u(t)  and  z(t) ,  we  find  that  the 
function  y(t)  can  be  expressed  as  a  function  of  r  =  v(t  -  tD)  by 
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y(t)  = 


-  5oI 

e  2V  cos  t 


nH 


n 


yQ 


n 


(2  +  n)v^o 


n  V 


T  2  V 


cos  t !  dr 


(22) 


where  yD  =  y(t  =  tQ)  and  tc  are  arbitrary  constants,  and  y(t)  is  a  general 
solution  to  the  following  equation: 


d  y  ,  /fj  ,  H  n\dy 
— 5  +  (H0  +  Hny  ;dt  + 

dt 


Mq  -  — r-  +  yn  + 

0  4  2  +  n 


H. 


(2  +  n) 


n  2  y2n 


y  =  0 


(23) 


where  n  is  an  integer. 

Solution  (22)  to  equation  (23)  was  first  given  by  Smith  (ref.  10), 
although  by  an  entirely  different  derivation.  Subsequently,  we  shall  use  this 
exact  solution  as  a  basis  of  an  approximate  solution  that  does  not  restrict 
the  nonlinear  restoring  terms  Mn  and  M2Q  to  the  values  given  by  equation  (21). 

It  is  interesting  to  see  the  explicit  forms  that  solution  (22)  takes  for 
specific  values  of  n.  Here  we  present  the  solutions  for  n  =  2  and  n  —  4, 
which  are  most  likely  to  be  of  practical  interest. 


n  =  2 


_  v 

e  2V  cos  t 


y  = 


H, 


HqT 


B  -  -r— —  e  V  [1  -  sin  9  sin(2T  -  cp)  ] 
4Hq 


(24) 


H, 


v2  =  MQ  -  -JJ—  ,  T  =  v(t  -  t0) 


sin  9  = 


H, 


n  =  4: 


y  = 


Hqt 

2.  V 


COS  T 


B  - 


H 

8a 


2H0t 


1  -  —  sin  \|r  sin(2T  -  ty) 


^  sin  cp  sin(4r  -  cp) 


(25) 


H. 


H, 


sin  = 


*  »< 


,  sin  cp  = 
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where  B  is  a  constant.  Formulas  (24)  and  (25)  are  useful  Because  they  indi¬ 
cate  explicitly  and  analytically  how  the  motion  may  be  affected  by  arbitrarily 
large  damping  nonlinearities. 

The  special  exact  solution  (22)  has  the  shortcoming  that  the  nonlinear 
static-moment  terms  ^  are  not  open  constants,  but  are  determined  by  the 
damping-moment  terms  H0  and  Hn  as  given  by  equation  (2l) .  We  would  like  to 
lift  this  restriction  and  find  approximate  solutions  for  arbitrary  values  of 
Mn-  To  do  this,  we  must  take  into  account  that  the  frequency  may  be  variable 
and  recast  the  equations  into  more  appropriate  forms. 
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dT 

dt 


(28c) 


=  v  1  +  — i-  f(a  cos  t ,  -av  sin  t)cos  t 
av2 

The  amplitude  variable,  a,  is  determined  as  a  function  of  t,  the  frequency 
variable,  by  solving  equation  (28b).  Using  this  solution,  a  =  a(r),  one  may 
then  find  the  relationship  between  t  and  the  original  independent  variable 
t  by  solving  equation  (28c)  in  the  form  of  a  quadrature 


(28d) 


Consider  now  the  application  of  the  modified  Kryloff-Bogoliuboff  equations 
to  our  particular  problem*  To  account  for  oscillations  for  which  _M0  <  0,  we 
shall  select  special  values  for  the  arbitrary  parameter  v2.  If  Mq  £  0,  then 
we  shall  choose  v2  =  |Mo| .  If  %  =  0,  then  we  shall  choose  v2  =  |Mnym  l> 
where  ym  is  a  suitable  characteristic  amplitude  of  the  motion.  Noting  that 
z  =  a  cos  t,  we  may  identify  the  function  f(a  cos  t,  -av  sin  t)  from  equa¬ 
tion  (26a)  as  follows: 


Expression  (29),  in  this  case,  is  not  in  convenient  form  to  solve  equa¬ 
tion  (28b)  because  the  variable  u  is  present,  and  u  must  be  found  as  a 
function  of  a(r)  and  t  by  combining  equations  (26b)  and  (28c)  to  eliminate 
the  variable  t.  Because  this  step  cannot  be  easily  accomplished,  we  instead 
introduce  a  new  variable ,  Q,  defined  as  follows: 

HqT 

Q  =  a^2  e  V  (30a) 

It  now  can  be  established  that  y  is  determined  as  a  function  of  t  by 

HqT 

y2(r)  =  Q(-r)e"  v  cos2  t  (30b) 

An  equation  for  Q(t)  now  remains  to  be  found.  From  equation  (30a)  it 
follows  that 

Hot  u 

,,2  da2  _  dQ  _  Q  dug  dt  _  gog  (30c) 

dT  dr  u2  dt  dr  v 

Substituting  equations  (26b),  (28b),  and  (28c)  into  equation  (30c),  we  obtain 
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_1_  dQ 
Q2  dT 


Hs 

2  +  s 


S-2 


sHqt 


H, 


2V  coss  t  +  ( ~  cos  t  +  2  sin  t 


f 

1  +  -  COS  T 

av2 


f 

av2Q 


(31a) 


This  equation  is  the  more  general  counterpart  of  equation  (19)  •  The  appropri¬ 
ate  expression  for  f  in  terms  of  Q  and  t  is  found  from  equation  (29)  to  he 


f 


av 


2 


+ 


(31b) 


Thus,  when  equation  (31b)  is  used  for  f  in  equation  (31a)  >  a  first-order  non¬ 
linear  equation  for  Q(t)  results. 

It  is,  in  general,  too  difficult  to  integrate  equation  (31a)  for 
Q  =  Q(t).  Nevertheless,  it  is  convenient  to  go  a  step  further  and  find  an 
integral  equation  for  y  that  will  serve  as  a  basis  for  approximate  solutions. 
We  first  note  that  f  can  be  expressed  in  terms  of  y2  by  using 
equation  (30b): 


Equation  (31a)  may  now  be  expressed  as 


(32a) 


(32b) 


1  IQ 

Q2  dT 


2  2 

—  COS  T 
V 


1  + 


+ 


COS 


T 


(33) 
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where 


r(t)  =  cos4  t  +  2  sin  t  cos3  t  (3*0 


Direct  integration  of  equation  (33)  gives 


(35) 


Substitution  of  equation  (35)  into  equation  (30b)  now  yields  an  integral  equa¬ 
tion  for  y(r): 


(36) 


Here  we  have  identified  the  constant  of  integration  B  as  l/y02>  where 

yD  =  y(T  =  0) . 

If  it  is  assumed  that  y  =  y(*r)  can  be  found,  then,  at  least  formally,  the 
corresponding  relation  for  t(r)  can  be  determined  from  equation  (28d): 


The  integral  equation  (36)  is  interesting  because  it  reflects  roughly  the 
nature  of  a  solution  that  one  might  expect  for  a  nonlinear  oscillation. 
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Furthermore,  when  v2  =  Mo>0,  Ms  =  0  for  s  >  1,  and  Hs  =  0  for  s  ^  0,2,  the 
exact  solution  (22)  for  n  =  2  is  recovered  immediately.  Equation  (36)  will 
thus  he  useful  for  a  perturbation  of  this  exact  solution. 


Equivalent  First-Order  Equation 


We  shall  now  derive  the  most  useful  and  most  important  equation  of  this 
analysis.  Since  y  =  uz,  the  first  derivative  of  y  can  be  written  in  terms 
of  u  and  z  as 


-  v/l  —  i  §2i\ 

dt  y\z  dt  u  dty 


(38) 


We  now  again  introduce  the  modified  Kryloff-Bogoliuboff  equations,  derived  in 
appendix  A,  so  that  z  =  a  cos  t  and  dz/dt  =  -av  sin  t*  In  addition,  we  make 
use  of  equations  (26b),  (28c),  and  (32b),  and  write  equation  (38)  as  follows: 


dy 

dT 


-y{  tan  t  +  i 


H, 


2  +  s 


s>o 


sin2  t  +  cos2  t 


M 


“S  1 

-  y 


s>o 


(39) 


The  general  second-order  differential  equation  (l)  with  t  as  the  inde- 
pendent  variable  has  now  been  reduced  to  a  single  first-order  equation  with  t 
as  the  independent  variable .  The  relation  between  t  and  r  is  again  provided 
by  equation  (37)*  Although  equation  (39)  is  equivalent  to  the  integral  equa¬ 
tion  (36)  because  it  describes  the  same  nonlinear  oscillation,  it  has  a  very 
different  appearance . 


Distortion  Due  to  Nonlinear  Damping 

Equation  (39)  can  immediately  be  used  to  determine  where  the  maximum 
amplitudes  occur  in  a  nonlinear  oscillation.  When  there  is  no  damping  in  an 
oscillation,  the  motion  is  periodic,  and  it  is  easy  to  deduce  that  the  maximum 
amplitudes  occur  at  Tm  =  nit,  where  n  is  an  integer.  For  a  general  oscilla¬ 
tion,  the  maximum  amplitude  occurs  when  dy/dr  vanishes,  and  thus  from  equa¬ 
tion  (39)  we  deduce  that 

tan  Tm  =  "  7  Yj  2_+— s  ymS 

s>o 


where  ym  =  y(%)  is  the  maximum  amplitude. 
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Formula  ( 40 )  indicates  how  much  the  maximum  amplitude  is  shifted  from  the 
point  of  symmetry,  Tm  =  nit.  Consider  a  special  simple  example,  namely,  the 
limit  cycle  for  the  exact  solution  (24).  The  limit  cycle  occurs  when  HQ  <  0, 
H2  >  0,  and  t  00  and,  in  this  case,  equation  (24)  yields 


The  value  of  ym  from  equation  (4l)  is 


and  equation  (4o)  now  yields  for  this  example 


and  the  limit  motion  is  distorted  from 
the  symmetry  position  by  an  amount 
cp  =  tan"1(H0/2v)  .  This  motion  is  shown 
in  figure  1  for  values  of  cp  =  0°, 

-15°,  -45°,  and  -75°-  When  |cp|  is 
large,  a  substantial  distortion  of  the 
nearly  cosine  oscillation  characteris¬ 
tic  of  small  nonlinearities  is  found, 
and  the  oscillation  is  no  longer  sym¬ 
metric  about  Tm  =  nrt. 


Figure  1.-  Distortion  of  motion  due  to 
nonlinear  damping. 


EXACT  SOLUTIONS  FOR  ZERO  DAMPING 


It  is  interesting  and  instructive  to  see  how  exact  solutions  for  zero 
damping  appear  in  the  present  formulation  of  nonlinear  oscillations.  Moreover, 
these  exact  solutions  will  be  useful  later  in  approximations  for  nonzero 
damping . 


General  Solution 


Let  us  now  consider  only  symmetric  motions,  so  that  s  is  even,  and  all 
the  damping  terms  are  zero  (Hs  =  0) .  The  motion  will  be  periodic  and  have  a 
maximum  amplitude  that  we  denote  by  ym*  We  introduce  a  normalized  form  of  y 
by  the  transformation 


In  addition,  we  introduce  the  coefficients  ms  defined  by: 

M  „  s 

(42) 

Note  that  if  MQ  /  0;  ms  is  a  ratio  of  the  magnitude  of  the  sth  nonlinear 
term  of  the  static  moment  to  the  magnitude  of  the  linear  term.  With  this  nota¬ 
tion,  one  may  express  equation  (39)  for  zero  damping  as  follows: 


rrn  = 


s  >  0 


dY0  _  -2Y0  tan 

d  t 

sin2  t  +  cos2  t 

s>° 

Equation  (43)  can  he  immediately  integrated  if  a  transformation  to  a  new 
independent  variable  x  =  tan  t  is  made.  If  the  boundary  condition  is  imposed 
that  Y0  =  1  when  t  =  nrt  (n  is  an  integer),  then 


(43) 


msYo 


S/2 


tan2 


2  +  s 


s>o 


(44) 


The  zero-damping  solution  becomes  complete  when  a  relationship  between  t 
and  t  is  established.  This  is  done  by  means  of  equation  (37) >  which  now 
becomes 
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1 

V 


T 


(45) 


t 


dT 


sin2  r1  +  cos2  t*  ^  msYo 
s>o 


s/2 


Formula  (45)  is  not  always  a  convenient  expression  to  use  because  YQ  is  gen¬ 
erally  an  implicit  function  of  t*  Consequently ,  an  alternate  expression  in 
which  t  is  determined  as  a  function  of  YQ  is  often  more  useful.  The  alter¬ 
nate  formula  may  be  found  by  changing  variables  according  to  equations  (43)  and 
(44) ,  and  thus  obtaining 


for 


(46a) 


nfl  <  t 


<  nit  + 


jt 

2 


and  0  <  YQ  <  1 


Expression  (46a)  has  been  obtained  by  more  conventional  methods . .  If 
m  =0  when  s  >  6,  then  equation  (46a)  can  be  evaluated  by  elliptic  integrals 
(see  refs.  11  and  12).  The  angular  frequency  w  of  the  oscillation  is  defined 

by 


(46b) 


An  approximate  formula  for  the  frequency  with  arbitrary  values  of  ms  has  been 
derived  in  reference  13  * 


Special  Case  of  a  Cubic  Static  Moment 


For  the 
easy  to  find 


case  of  a  linear -plus -cubic  moment,  ms  =  0  when  s  >  4,  it  is 
an  explicit  function  for  Yq(t)  from  equation  (44);  namely, 


IT 


Yo  - 


(n/m2)cos2  t 


sin2 


t  +  m0  cos2  t  +  ^sin2  t  +  m0  cos2  t)  +  |i  cos4 


(47) 


where  n  =  m2(2m0  +  m2)  . 

With  this  explicit  representation,  equation  (37)  yields  the  following 
relation  between  t  and  t 


t  '  t0  =  T 


dT' 


o 


J( sin2  t*  +  mQ  cos2  tt)2  +  P-  cos 


(48a) 


4  t* 


It  follows  that  the  frequency  is  given  by 


v 

co 


Jt/2 


dT 


J( sin2  t  +  mQ  cos2  t)2  +  [x  cos4  t 


(48b) 


It  is  interesting  to  compare  equation  (48b)  for  the  frequency  with  the 
corresponding  formula  (46b) ,  which  may  be  evaluated  in  terms  of  elliptic  inte¬ 
grals.  The  quadrature  (48b)  is  a  proper  integral  and  may  be  evaluated 
numerically  without  difficulty. 


Formulas  (48)  contain  four  possible  cases  of  interest,  depending  on  the 
value  of  the  moment  terms  Mq  and  M2.  To  discuss  them,  we  adopt  the  following 
convention.  If  Mn  is  positive,  it  is  called  stable;  if  it  is  negative,  it 
is  called  unstable.  _The  various  possibilities  are  stated  by  describing  the 
term  K 0  and  then  M2.  For  a  cubic  static  moment  we  thus  have  the  following 
cases :2 


Stable- stable  and  stable-unstable  static  moment 


Mq  >  0,  v2  =  Mq,  m0  =  1 


0  <  m2  <  oo 

-1  <  m2  <  0 


( stable- stable ) 

( stab  le  -un  stab  le ) 


Unstable- stable  static  moment 
Mo  <  0,  v2  =  |Mq|,  mQ  =  -1 
Pure  cubic  static  moment 

Mq  =  o,  V2  =  M2ym2  >  0,  mQ  =  0,  m2  =  1 


2  <  m2  <  oo 


>  (49) 


2Note  that  if  m2  <  -1  for  the  stable -unstable  case ,  tumbling  motion 
results;  also,  if  1  <  m2  <  2  for  the  unstable -stable  case,  unsymmetric 
motion  about  a  trim  angle  occurs.  Since  neither  of  these  cases  is  covered  in 
the  present  analysis,  initial  conditions  must  be  such  that  the  restrictions  on 
m2  are  .observed. 
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The  oscillations  (eq.  (4-7))  for  the  cubic  moments  are  plotted  in  figures  2(a), 
(b),  and  (c)  to  show  the  effect  of  m2  for  the  various  cases.  In  these  fig¬ 
ures,  =  y/ym  is  plotted  versus  the  transformed  independent  variable,  t, 

over  1  quarter  cycle  of  oscillation.  As  indicated  by  equation  (42),  m2  is  a 
ratio  of  the  magnitude  of  the  nonlinear  part  of  the  static  moment  to  the  magni¬ 
tude  of  the  linear  part.  The  value  m2  =  0  represents  the  pure  linear  case, . 
m2  =  oo  the  pure  cubic  case.  In  figure  2(d) ,  y/ym  plotted  versus  the  orig¬ 

inal  independent  variable ,  t,  normalized  in  such  a  way  that  the  curves  have  the 
same  frequency.  This  is  equivalent  to  comparing  each  nonlinear  moment  with  a 
different  linear  moment,  in  each  case  the  linear  moment  being  chosen  to  give 
the  same  frequency  of  oscillation  as  did  the  particular  nonlinear  moment.  With 


(c)  Variation  with  t;  unstable -stable  1-3  (d)  Variation  with  t. 

static  moment. 

Figure  2.-  Effect  of  m2  on  oscillatory 
behavior • 


this  kind  of  plot,  it  can  be  seen  that  most  nonlinearities  in  the  moment  have 
no  drastic  effect  on  the  appearance  of  the  motion  (although  the  frequency  is 
greatly  affected).  For  instance,  all  stable-stable  cases  would  fall  within  the 
shaded  region  of  figure  2(d).  Furthermore,  the  vast  majority  of  stable- 
unstable  cases  ever  encountered  would  be  bracketed  by  the  curves  labeled  0.0 
and  -0.8;  most  unstable -stable  cases  would  be  bracketed  by  the  curves  labeled 
2.5  and  00 . 


APPROXIMATE  SOLUTIONS 


We  now  wish  to  find  solutions  for  arbitrary  values  of  Hs  and  Ms.  With 
the  use  of  the  special  exact  solutions  already  discussed,  it  is  possible  to 
find  approximate  solutions  valid  when  the  parameters  Hs  and  Ms  are  small  in 
a  certain  sense.  Essentially,  we  shall  find  first-order  perturbation  solu¬ 
tions  to  the  exact  solutions  that  are  special  cases  of  equations  (36)  and  (39) • 
We  shall  consider  only  symmetric  motion  and  hence  make  the  restriction 
Hs  =  Ms  =  0  when  s  is  odd. 


Perturbation  Solution  Using  the  Integral  Equation 

Because  the  integral  equation  (36)  has  an  exact  solution  for  a  special 
case,  that  is,  when  Mq  >  0  and  Ms=0,  s>0,  and  when  Hs  =0,  s  >  4,  it  is 
natural  to  ask  if  there  exists  a  small  perturbation  of  this  exact  solution  in 
which  certain  parameters  are  small  in  an  appropriate  sense •  The  exact  solu¬ 
tion  can  be  expressed  in  the  form  of  equation  (22),  with  n  =  2, 


y2  ( T) 
(o) 


.  2oi 

e  v  cos2  t 


dT1 


(50) 


Expression  (50)  can  he  considered  as  a  zeroth  order  approximation.  Since  equa¬ 
tion  (36)  is  an  integral  equation,  it  is  straightforward  to  use  equation  (50) 
and  determine  a  first-order  approximation  to  equation  (36)  as  a  first  step  in  a 
scheme  of  successive  approximations. 


Consistent  with  the  zeroth  approximation,  we  require  that 

that 


M0  >  0,  and 


Ms  vs 
v2  7 


S>2 


«  1 


(51a) 


and 


I 

S>4 


^Hs 

(2  +  s)H2 


S-2 

y 


«  1 


(51b) 
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With  these  conditions,  we  can  expand  the  denominator  of  the  integrand  of  the 
integral  in  equation  (36)  in  a  Taylor  series  and  keep  the  first  term.  The 
first-order  approximation  can  now  be  obtained  by  replacing  the  ^y2  that 
appears  under  the  integral  by  the  zeroth-order  approximation  y^Q)  retain¬ 

ing  terms  consistent  with  equations  (51) .  The  first-order  approximation  thus 
can  be  written  as 


e 


2 

cos 


T 


H0t' 


S>2 


Hs  S-2 

2^y(o) 


S>2 


M 


S-2 

(o) 


f 


R(t') 


Hg 

2v 


y(o) 


cos 


d-r* 


(52) 


Formula  (52)  for  the  first  approximation  is  expressed  in  terms  of  quadra¬ 
tures  involving  the  zeroth-order  approximation.  In  general*  these  quadratures 
are  difficult  to  evaluate  in  terms  of  tabulated  functions.  For  the  special 
case  when  the  nonlinear  damping  terms  are  zero*  however*  the  quadratures  can  be 
evaluated  with  the  use  of  elementary  functions.  We  now  write  the  corresponding 
formal  solution  for  the  frequency  function  t  as  a  function  of  t.  Starting 
with  equation  (37)  and  using  the  approximations  consistent  with  the  first- 
order  solution*  we  obtain 


v(t  -  t0) 


(53a) 


where*  as  before*  y^  is  given  by  equation  (50)  •  Since  higher  order  products 

of  M  are  neglected  in  the  first  approximation*  we  can  also  express  t 
explicitly  in  terms  of  t  by  rearranging  equation  (53a)  to  read 


T  —  v(t  ”  tQ) 


(53b) 


In  many  practical  problems  of  inheres!  in  aerodynamics,  the  nonlinear 
restoring-moment  terms  are  not  small,  and  thus  the  above  first-order  approxima¬ 
tions  are  not  particularly  useful*  The  damping,  however,  is  usually  small; 
consequently,  results  having  wide  application  can  be  obtained  from  a  perturba- 
tion  solution  of  the  first-order  differential  equation  (39) . 
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Perturbation  Solution  Using  the  First-Order 
Differential  Equation 


Let  us  write  equation  (39)  in  a  normalized  form  that  naturally  suggests  a 
perturbation  procedure.  We  introduce  a  characteristic  amplitude  ya  that  we 
shall  identify  later  and  define  the  following  symbols: 


(54) 

s  even 


Since  Ms  &  Ms,  ms  is  essentially  a  ratio  of  the  magnitude  of  the  sth  non¬ 
linear  term  of  the  static  moment  to  the  magnitude  of  the  linear  term;  hg  is  a 
similar  ratio  for  the  damping  moment  but  multiplied  by  a  constant,  l/(2  +  s) . 
With  this  notation,  equation  (39)  becomes 


Since  an  exact  solution  exists  for  equation  (55)  when  h  =  0,  namely, 
equation  (44),  we  seek  a  perturbation  solution  of  equation  (55)  valid  when 
h  «  1.  Therefore,  we  consider  that  h  «  1  and  hs  =  0(l)  and  assume  over  a 
half  cycle,  at  least,  that  Y(t)  can  be  expanded  in  a  series  of  the  following 
form: 


00 


Y(t) 


-I 


hmYm(x) 


m=o 


(56) 


Substituting  equation  (56)  into  (55)  >  collecting  terms,  and  requiring  that  the 
coefficient  of  each  power  of  h  vanish  identically  yields  an  infinite  set  of 
equations  for  the  functions  Ym-  We  shall  confine  our  interest  to  the  first- 
order  perturbation  and  thus  to  only  the  functions  Y0  and  Y1.  The  equations 
for  Y0  and  Yx  are  found  to  be 


dYQ  2Y0  tan  t 

d7  =  jTq 


(57a) 
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D 


dYi  +  2  tan  t 


°  dr  D 


sin2  t  +  cos2  t  )  (  1  -  )  m„Y. 


s>o 


S/2 
S~0 


Y  =  -2Y  )  h  YS/2 

o  /  s^o 


S>0 


(5Tb) 


where 


S  /  2 

Do  =  sin2  t  +  cos2  t  )  msY0 

s>o 


(57c) 


With  the  use  of  equations  (57a)  and  (57c),  it  is  not  difficult  to  verify  that 
an  integrating  factor  for  equation  (57b)  is  merely  sec2  t,  so  that  equa¬ 
tion  (57b)  can  be  e:xpressed  as 


sec 


Yx)  =  -2  sec2  t  Y0  ^ 


,S/2 

O 


s>o 


(58) 


With  the  above  equations  we  will  be  interested  in  the  solution  between  two 
consecutive  maximum  amplitudes,  say  the  nth  maximum  amplitude,  y  ,  and  the 
(n+i)th,  yn+1.  The  zeroth  order  equation  (57a)  yields  a  periodic  solution 
with  a  given  maximum  amplitude .  We  now  select  ya  so  that  the  maximum  ampli¬ 
tude  squared  of  the  zeroth  order  solution  is  the  mean  of  the  squares  of  the 
nth  and  (n+i)th  maximum  amplitudes 


y  2  =  i(y 

2  ■ 


2  +  y|+i) 


We  can  now  take  equation  (44)  to  be  the 
solution  of  (57a)  and  identify  ym 
with  ya.  The  situation  is  shown  in 

sketch  (a) . 

Correct  within  the  first-order 
approximation,  yn  occurs  at  t  =  ntf 
and  yn+1  occurs  at  t  =  (n  +  l)rt* 

We  integrate  equation  (58)  so  that 
Y^mt)  =  -Yx(nJt  +  it) 


Yi(t)  =  - 


2  COS2  T 


sec 


(n+|)it 


I 

s>o 


*sYo 


S/2 


dT* 


(59) 
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The  function  Y1(t)  is  expressed  by  a  quadrature.  Since  Yd(t)  can  he  written 
explicitly  as  a  function  of  t  only  for  special  cases,  such  as  the  cubic 
restoring  moment  (eq.  (47)),  a  general  explicit  form  for  the  integrand  cannot 
be  found.  With  a  change  in  variable,  however,  the  quadrature  can  be  expressed 
in  the  form  of  elliptic  or  hyperelliptic  integrals . 

Consider  a  change  of  variable  in  the  integral  of  equation  (59)  to  Y0,  as 
defined  by  equations  (57a)  and  (44) .  Then,  in  terms  of  YQ,  the  function  Y1 
can  be  written  as 


Yx(t) 


-Y0G(Y0) 


(60a) 


for  nit  <  t  <  (n  +  l)fl.  The  plus  sign  is  used  when  [n  +  (l/2)]it  <  t  <  (n  +  l)it 
and  the  minus  sign  when  nit  <  t  <  [n  +  (l/2)]it.  It  is  clear  that  Y^t)  and 

G(Y0)  are  odd  functions  of  [r  -  [n  +  (l/2)]it}.  When  the  polynomial  in  the 
radical  of  equation  (60b)  is  a  quartic  or  less,  then  the  integral  can  be  eval¬ 
uated  in  terms  of  elliptic  integrals  and  elliptic  functions.  We  shall  eval¬ 
uate  explicitly  the  complete  integral  later  for  special  cases  of  the  damping 
decrement . 

The  corresponding  first-order  form  of  the  frequency  equation  (37)  can  now 
be  determined.  Inserting  equation  (56)  into  (37)  and  keeping  only  first-order 
terms  in  h  yields 


where  here  t  =  t0  when  t  =  [n  +  (1/2) 3 « .  This  completes  the  formal  first- 
order  solution.  Because  the  first-order  integral  in  equation  (6l)  is  very  com¬ 
plicated,  we  shall  not  investigate  it  further  here,  but  merely  note  that  it  can 
be  integrated  numerically  without  great  difficulty  when  Y0  and  Y1  are  known 
functions  of  t. 


THE  DAMPING  DECREMENT 


We  now  wish  to  derive  formulas  that  portray  the  distinguishing  features  of 
the  nonlinear  oscillation.  Of  an  oscillation  that  is  sufficiently  well-defined, 
two  distinguishing  features  are  the  maximum  amplitudes  and  the  half  period. 

(By  half  period,  we  mean  the  time  interval  between  two  successive  maximum  ampli¬ 
tudes  of  the  motion.  In  this  sense,  it  is  a  pseudo-period,  since  the  motion  is 
not  necessarily  periodic.  Given  the  half  period,  we  can  likewise  define  an 
angular  frequency.)  A  pertinent  parameter  to  consider  is  the  ratio  of  two  suc¬ 
cessive  maximum  amplitudes,  called  the  damping  decrement.  A  slightly  more  sig¬ 
nificant  parameter  is  the  logarithmic  decrement  because  it  is  a  constant  for  a 
strictly  linear  oscillation  and  thus  forms  a  good  standard  for  comparison. 


Decrement  for  Small  Damping 

Consider  now  the  first-order  solution  that  neglects  terms  of  order  h  in 
equation  (56) j  that  is, 

Y(t)  =  Yd(t)  +  hY  x( t)  +  0(h2)  (62a) 


and  yn+1  to  be  given  by 

1  +  hY^nit)  +  0(h2)  (62b) 

1  +  hYx(nrt  +  Jt)  +  0(h2)  (62c) 

Noting  that  Y1(nrt)  =  -Y1(nit  +  it),  dividing  one  equation  by  the  other,  and 
expanding  for  small  h,  we  find  the  damping  decrement  to  be 

=  1  -  2hY1(nit  +  it)  +  0(h2)  (63) 

Taking  the  logarithm  of  this  expression  and  again  expanding  for  small  h 
yields  the  following  formula  for  the  logarithmic  decrement: 


From  this  expression  we  find  yn 


In 


Yn 

Yn+i 


=  -hY1(nit  +  it)  +  0(h2) 


(64) 


25 


Here  Y1(rut  +  Jt)  can  be  evaluated  by  equation  (59)  or  (60) .  Before  proceeding, 
however,  it  is  convenient  to  modify  equation  (64)  by  defining  an  effective 
linear  damping. 


Effective  Linear  Damping 

In  a  strictly  linear  oscillation  it  is  known  that  the  logarithmic  decre¬ 
ment  is  a  constant  with  the  value 

H&l-^  (65 

where  v  is  the  frequency  defined  hy  v2  =  | M0 1 .  Analogously ,  we  can  define 
an  effective  linear  damping  HQ  by 


yn+i 


H0  it  Hn 

°e  _  ^  jt  v  °e 

"S3"  "  n  2  <3 


where  now,  consistent  with  the  first-order  analysis,  we  define  the  frequency 
w  by 


v  =  1 

0)  Jt 


sin2  t  +  cos2  t  )  msYo/2 


(67a) 


v  =  1 

CjO  Jt 


Yo  7,  rfs  1 '  Y° 


(67b) 


^he  logarithmic  decrement  can  be  written  in  more  conventional  form  if  we 
introduce  the  so-called  damping  ratio  &  as  6  =  H0/2MQ.  Then,  realizing  that 

v  =  =  J%  -  (H02/4)  for  a  linear  oscillator ,  we  have 


Tn+1  I  s/1  -  & 

For  a  linear  system,  the  logarithmic  decrement  is  independent  of  amplitude, 
and,  therefore,  we  can  also  write 


1  - 
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We  now  combine  equations  (64)  and  (66)  to  get 


^  =  _  2  co  Y  (nJt  +  „)  +  0(h)  (68) 

H0  «  v  i 

Expression  (68)  is  an  important  formula,  and  we  shall  investigate  it  in  some 
detail . 


General  Formula  for  Effective  Linear  Damping 


The  formal  first-order  value  for  the  effective  linear  damping  can  he 
found  by  evaluating  Y  (nit  +  it)  in  formula  (68)  by  (60)  (noting  that 


Y0(nit  +  it)  =  l) 


H, 


2  “ 
ft  V 


H0 


m, 


s>o 


j 


A 


m. 


S+2 


2  +  s 


2  +  sY 


o 


s>o 


■) 

/s>o 


S/2 

hsY0  dY( 


(69) 


where  co/v  is  given  by  equation  (67b) •  Formula  (69)  can  be  explicitly  evalu¬ 
ated  in  terms  of  tabulated  functions  for  restricted  combinations  of  ms.  We 
shall  investigate  these  special  solutions  in  order  to  gain  an  appreciation  for 
the  influence  of  the  nonlinearities  on  the  damping,  and  to  ascertain  the  asymp¬ 
totic  values  of  the  general  formula  in  various  limits.  On  the  basis  of  this 
information,  an  approximate,  but  quite  general,  formula  can  be  established 
that  is  valid  for  almost  the  whole  range  of  possible  values  of  the  nonlinear 
parameters . 


Cubic  Damping  and  Cubic  Static  Moment 


Formula  (69)  can  be  evaluated  in  terms  of  elliptic  integrals  for  a  cubic 
static  moment;  that  is,  ms  =  0  when  s  >  4.  For  the  sake  of  simplicity,  we 
also  consider  only  a  cubic  damping  moment;  that  is,  Hs  =  0  when  s  >  4.  For 
this  case,  formula  (69)  appears  as 


H 


o 


1  W 
ft  V 


m0  +  m2 


m  p 

HU  h — 2~ 


1  +  Y, 


■i 


(1  +  2h2Y0)  dY, 


Yn( 1  -  Y0) 


mr 


+  *#(1  +  Yft) 


(70a) 


where 
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1 


V 

CO 


dY, 


fY0(l  -  Y0) 


m2  /  X 

mQ  +  -f  (1  +  Y0) 


(70b) 


Two  representations  are  required  for  the  various  static  moments  defined 
ty  equation  (49)  to  express  (70a)  in  terms  of  elliptic  integrals.  (No  such 
distinction  is  needed  for  the  various  cubic  damping  moments  that  could  exist.) 
These  are  as  follows: 


H0e  2mQ  +  m2  H  f  m2  7£\  (y±  m2  7J\ 

Hq  mQ  +  m2  [_2  \  +  2m0  +  m2  Ky  2  VK  ^o  +  m2  K/. 


(71a) 


K  =  K(k) 


Jt/2 


dcp 


Jl  -  k2  sin2  cp 


E  =  E(k) 


sin2  cp  dcp 


Complete  elliptic  integral 
of  first  kind 

(71b) 

Complete  elliptic  integral 
of  second  kind 


(7ic) 


Stable-stable  and  unstable- stable  static  moment:  (formula  236. 16,  ref.  12) 


7,  s  h  [(k2  -  1)K  +  E] 

1  k 

7  =— 4  t(3k4  -  5k2  +  2)K  +  2 (2k2  -  l)E] 

2  3k 

=  -t?  [4(2k2  -  1)7  +  3(1  -  k2)7  ] 

3  5k  2 

2  1  m2 

K  2  mQ  +  m2 


(72) 


Stable-unstable  static  moment: 


(formula  23^*16,  ref.  12) 


7X  =  h  (K  -  E) 
x  k 

7  =  “^4  [(2  +  k2)K  -  2(1  +  k2)E] 

2  3k4 

7 3  =  ^(k2  +  1)72  -  37  J 

5k 

k2  = 

2  +  m2 


(73) 


An  equivalent  expression  for  HQe/H0  for  a  cubic  static  moment,  but  a 
linear  damping  moment,  that  is,  h2  =  0,  has  been  given  by  Murphy  (ref.  5)  and 
Murphy  and  Hodes  (ref.  6).  Murphy  obtained  his  results  in  a  different  manner. 
He  assumed  that  the  damping  was  small  and  then  used  the  exact  elliptic  func¬ 
tions  obtained  for  zero  damping  in  an  averaging  process  over  a  cycle  of  motion. 
The  perturbation  method  of  Murphy  thus  yields  the  same  value  for  the  effective 
damping  as  the  present  analysis  for  a  cubic  static  moment,  but  it  is  difficult 
to  extend  the  method  of  Murphy  to  higher  order  static  moments  owing  to  the 
complexities  in  representing  more  complicated  static-moment  solutions  by  ellip¬ 
tic  functions.  An  approximate  formula  for  Hoe/H0  for  a  cubic  static  moment 
and  linear  damping  has  also  been  derived  by  Murphy  (ref.  7)  by  a  "quasi-linear" 
method,  and  we  will  discuss  this  later  in  the  section  entitled,  "Comparison  for 
a  Cubic  Moment." 


Expression  (71a)  for  Hoe/H0  in  terms  of  elliptic  integrals  is  actually 
a  fairly  involved  formula.  Also,  the  quadratures  given  by  equation  (70)  are 
difficult  to  evaluate  numerically  because  of  their  singular  nature  at  both 
limits.  It  is  interesting  that  the  effective  linear  damping  can  be  expressed 
by  another  set  of  quadratures  that  possibly  may  be  easier  to  evaluate  numeri¬ 
cally  (especially  with  electronic  computers)  than  either  equation  (70a)  or 
(71a) .  We  obtain  this  alternate  formula  by  transforming  equation  (70a)  from 
the  Y0  variables  to  the  r  variables  with  the  use  of  equations  (43),  (44), 
and  (47)  .  The  alternate  result  is 
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H, 


Qg  _  2w  gmo+ma 
H  Jtv  m0+m2 


+  2h2(2m0+m2) 


•Jt/2 


dT 


(sin2  T+m0  cos2  T)+N/(sin2  T+m0  cos2  t)2+^  cos4  t 


■at/s 


COS  T  dT 


(sin2  T+m0  cos2  T)+N/(sin2  T+m0  cos2  t)2+M-  cos4  t 


(74a) 


where 


29 


dr 


v 

U) 


i  f 


Jt/2 


^(sin2  t  +  m0  cos2  t)2  +  P  cos4  t 


(74b) 


and 


|i  =  m2(2m0  +  m2) 


The  above  quadratures  are  proper  integrals  and  can  be  evaluated  numerically 
without  difficulty.  It  is  worth  mentioning  that  equation  (74a)  can  easily  be 
extended  to  higher  order  damping  terms. 

There  are  several  ways  in  which  formulas  (7la)  and  (74a)  can  be  plotted. 

If  one  is  interested  in  knowing  H0e/H0  for  given  values  of  h2  and  m2,  carpet 
plots  are  useful.  Such  carpet  plots  are  shown  in  figures  3(a),  (b),  and  (cj. 
Because  of  the  ease  of  making  linear  interpolations,  these  plots  are  convenient 

to  use. 


(a)  Stable -stable  1-3  static  moment. 


(b) 


Stable -unstable  1-3 


static  moment. 


(c)  Unstable -stable  1-3  static  moment. 


Figure  3»-  Carpet  plot  of  effective  linear 
damping. 
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Another  method  of  plotting  illustrates  how  /Hq  varies  with  amplitude 
for  given  values  of  the  damping-  and  static-moment  parameters.  We  thus  plot 
H0e/H0  as  a  function  of  m2  and  set  h2  =  (h2/m2)m2  so  that 

ha  _  _Hg  vf  =  l^ol 
m2  4hg  M2  4hq  m2 

can  he  varied  as  a  parameter  that  does  not  depend  on  amplitude  for  a  given 
oscillator.  These  plots  are  illustrated  in  figures  4(a),  (b),  and  (c).  The 
curves  denoted  by  h2/m2  =  0  were  given  by  Murphy  and  Hodes  (ref.  6).  The 
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(a)  Stable -stable  1-3  static  moment. 


(b)  Stable-unstable  1-3  static  moment. 


prig 


(c)  Unstable -stable  1-3  static  moment. 
Figure  4.-  Effective  linear  damping  versus  m2. 
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curves  in  figures  4(b)  and  (c)  are  singular  for  the  stable -unstable  moment  at 
xa2  s  M2y  2/Mq  =  -1  and  for  the  unstable-stable  moment  at  m2  -  M2yn  / | M0 |  =  2. 
The  dashed  lines  are  the  asymptotes  for  m2  00 •  It  is  possible  to  show  by  an 
expansion  for  large  m2  that  the  asymptotes  are  given  by 


where 


or 


mOh2 

m2 


m2 


T  =  gamma  function 


0.667  +  0.296 


mQh2 

m2 


+  0.7311 


(75) 


Likewise,  for  the  stable-stable  and  stable -unstable  moments,  the  variation  as 
m2  -*•  0  is  given  by 


Ho 


1  + 


m2  +  0(m22) 


(76) 


Another  interesting  aspect  of  these  curves  is  the  appearance  of  maxima  and 
minima  in  some  cases.  If  these  curves  were  established  by  experiment,  this 
behavior  might  be  considered  perplexing  without  theoretical  background.  They 
portray  the  effects  of  strong  interactions  of  the  nonlinear  damping  and  static 
moments.  Still  another  interesting  feature  is  the  curve  labeled  h2/m2  =  5/4 
in  figure  4(b).  Detailed  examination  of  equation  (74)  shows  that  this  value 
is  the  dividing  point  between  cases  where  a  limit  cycle  can  or  cannot  exist, 
given  a  stable-unstable  static  moment.  For  h2/m2  <  5/4,  there  is  no 
possibility  of  a  limit  cycle. 

All  curves  shown  in  figure  4  were  compared  with  results  obtained  from 
numerical  integrations  of  equation  (l)_^  In  these  numerical  integrations,  var¬ 
ious  values  of  the  parameter  h  s  H0/|M0|1/2  were  chosen  up  to  a  value  of 
about  0.1.  This  encompasses  the  practical  range  of  dairping  encountered  in 
ballistic-range  testing.  In  all  cases,  the  results  of  the  numerical  integra¬ 
tions  agreed  almost  exactly  with  the  results  from  formula  (71a)  or  (74a)  shown 
in  figure  4.  To  the  scale  these  curves  are  plotted,  no  differences  can  be 
shown . 

When  H0  <  0,  a  limit  cycle  exists  when  H0g  =  0.  The  amplitude  of  the 
limit  cycle  may  be  found  by  setting  H0g  =  0  in  equation  (7la)  or  (74a). 

Using  equation  (74a)  we  get 
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o  '  l 


Jt/2 
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(sin2  t+di0  cos2  T)+„/(sin2  T+m0  cos2  r)2+h  cos4  t 


2(2mr,+m2)  C^!2 


cos2  T  dT 


L  [( 


sin2  T+m0  cos2  t 


)+J( sin2  T+m0  cos2  t)2+M-  cos4  t_ 


(77) 


An  alternate  formula  is  obtained  from 
equation  (71a),  which  is  equivalent  to 
the  result  found  by  Murphy  and  Hodes 
(ref.  6)  by  investigating  the  singu¬ 
lar  nature  of  the  limit  cycle  in  an 
amplitude  plane*  The  effect  of  the 
cubic  static  moment  on  the  amplitude 
of  the  limit  cycle  is  demonstrated  in 
figure  5 •  Except  in  the  region  of  the 
singularities,  the  effect  of  m2  is 
not  large,  yielding  a  value  of 
(-h2)j  =  0.9118  as  m2  -*•  00  contrasted 
to  the  well-known  value  ( _h2)  7,  =  1 
when  m2  -*•  0.  It  will  be  shown  later 
that  for  quintic  and  higher  order 
static  moments,  the  effects  are  much 
greater . 


Arbitrary  One -Term  Static  Moment 


The  quadratures  in  equation  (69)  can  be  evaluated  in  terms  of  tabulated 
functions  for  the  special  case  of  an  arbitrary  one-term  static  moment .  This 
case  yields  the  asymptotic  values  that  are  obtained  as  the  highest  order  term 
in  an  arbitrary  polynomial  static  moment  tends  to  large  values. 


We  now  consider^  Mg  ^  0  for 
we  choose  v2  =  Mnyn  >  0  so  that 
now  appears  as 


s  =  n,  where  n  is  an  even  integer.  Hence 
ms  =  0,  s  n,  and  mn  =  1.  Expression  (69) 
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and 


v 

0) 


1 

jt 


(78b) 


( n+2 )  /  2. 

With  the  substitution  x  =  Y0  ,  these  quadratures  can  be  written  in  the 

form  of  beta  functions,  which  in  turn  can  be  evaluated  in  terms  of  gamma 
functions,  T.  For  equation  (78a),  we  thus  obtain 


where 
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s,  n  >  0 


(79b) 


Although  formula  (79)  was  derived  for  a  one-term  static  moment,  it  also 
represents  the  asymptotic  value  of  HQ  /H0  for  a  polynomial  static  moment  in 
which  the  nth  term,  mn,  is  dominant  and  stable  and  tends  to  infinity.  For 
purely  linear  damping  it  is  interesting  that  the  asymptotic  value  is  simply 
H°e/Ho  "  k/[\  +4),  On  the  other  hand,  for  a  linear  static  moment,  n  =  0, 
equation  (79a)  can  be  written 


H, 


H, 


1+2 


V  ±1?  • 

2*4*6 


(s  -  1) 
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S>2 


(79c) 


Redd  et  al.  (ref.  9),  considering  only  a  linear  static  moment,  computed 
the  same  value  as  equation  (79c)  up  to  s  =  8  by  a  different  method.  Their 
conclusion,  however,  that  a  nonlinear  static  moment  does  not  significantly 
alter  the  formula  (79c) >  except  near  unstable  trim,  is  not  substantiated  by  the 
present  analysis.  Here  it  is  found  that  equation  (79c)  is  a  special  case  of  a 
more  general  formula  (79a) ,  and  each  type  of  nonlinear  static  moment  may 
indeed  affect  H0e/H0  significantly. 

Expression  (79a)  may  be  used  to  examine  the  effects  of  the  nonlinear 
static  moment  on  the  value  of  the  limit  amplitude.  Consider  a  two-term  damp¬ 
ing  moment  given  by  Hs  =  0  when  s  /  0  or  k  (k  even)  so  that  the  first  term 
is  linear  and  the  second  term  is  nonlinear  of  arbitrary  order.  The  limit  cycle 
occurs  when  Hq  <  0  and  may  be  determined  by  setting  H0g  =  0  to  obtain 
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(80) 
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(n  +  4)7kn 


The  value  of 
values  of 


(-h^)^  is  a  function  of  both  k  and  n,  and  some  of  the 
are  given  in  the  following  table: 


\  i 
n  1 
0 

;  2 

4 

6 

8 

1.000 

1-333 

1.600 

1.829 

2 

.912 

1.167 

1.368 

1.540 

4 

.870 

1.084 

1.250 

1-392 

6 

.845 

1.034 

1.178 

1.300 

8 

.829 

1.001 

1.130 

1.238 

00 

-750 

.833 

-875 

.900 

It  can  be  seen  from  this  table  that  the  value  of  (“h^)}  decreases  for  a 
given  value  of  k  as  n  increases-  Here  we  recall  that  n  =  0  represents  a 
purely  linear  static  moment. 


Expansion  for  Small  Static-Moment  Nonlinearities 

Consider  the  value  of  H0e/H0  when  the  nonlinearities  in  the  static 
moment  are  small.  Here  we  assume  MQ  >  0  so  that  mQ  =  1,  and  we  assume  that 
ms  is  small  so  that  higher  order  products  and  cross-products  can  be  neglected. 
With  these  conditions,  it  is  possible  to  expand  equation  (69)  for  small  ms 
and  evaluate  each  integral  term  by  term.  The  result  can  be  written  as  follows 
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Some  of  the  values  of  b^s 
follows : 


for  various  values  of  k  and  s  are  tabulated  as 
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Expression  (76)  for  the  cubic  damping  and  cubic  static  moment  is  a  special 
case  of  (8la) .  Formula  (8la)  will  be  of  use  in  the  construction  of  a  general 
approximate  formula  for  H0e/HQ. 


APPROXIMATE  GENERAL  FORMULA.  FOR  H^/Hq 


Although  the  formal  representation  (69)  is  a  very  general  result,  it  is 
difficult  to  evaluate  for  arbitrary  combinations  of  the  damping  and  static 
moments.  Even  for  the  cubic  case,  the  evaluation  in  terms  of  elliptic  inte¬ 
grals  is  rather  complicated.  Based  on  the  special  cases  (79)  and  (8l),  how¬ 
ever,  it  is  possible  to  construct  an  approximate  formula  that  is  extremely 
accurate . 

A  perusal  of  the  formal  solution  (69)  indicates  that  H0e/H0  can  be 
expressed  by 

le  _  V 

Ho  L 

k>°  ^  pksms 

S>0 


V“ 


s  s 


S>0 


(82) 


where  aks  and  |3ks  are  functions  of  the  parameters  m0,m2,m4,.  .  .  .  The 
functions  aks  and  (3ks  depend  only  weakly  on  these  parameters,  however,  and 
thus  we  wish  to  find  approximate  values  of  akg  and  Pks  that  are  independent 
of  m0>m2,m4,.  ...  We  can  do  this  "by  requiring  that  equation  (82)  have  the 
proper  behavior  when  the  nonlinear  static  moment  is  very  large  and  stable  and 
when  the  nonlinear  static  moment  is  very  small;  that  is,  equation  (82)  should 
have  the  same  behavior  as  (79)  and  (8l)  in  these  particular  limits. 
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Formula  (82)  is  thus  guaranteed  to  be  valid  when  the  dominant  term 
m0jm2,m4,  .  .  .  00  and  also  when  the  terms  m0>m2,m4,  .  .  .  0. 


We  begin,  somewhat  arbitrarily,  by  requiring  that  f^o  =  1.  We  do  this 
in  order  to  agree  with  the  leading  term  of  the  sum  £  ms  that  appears  in  the 

denominator  of  equation  (69)  •  It  can  be  shown  that  if  equation  (82)  is  to 
reduce  to  both  (79)  and  (8l)  in  these  respective  limits,  then  a^s  and 
are  given  by 


<\Q  =  2l(k  -  1)  ,  Pko  =  1 


‘‘ks 


^ks^ks 
“ko  -  ?ks 


(83) 


ks 


kS  -  ?ks  ' 


s  >  2 


where  7ks,  bks,  and  l(n)  are  defined  by  equations  (79b)j  (8lb) ,  and  (8lc). 
Some  values  of  a^s  and  3ks  are  presented  in  the  following  tables.4 


^s 


0 

1 

2 

3 

4 

5 

6 

7 

8 

0 

2.0000 

1.2732 

1.0000 

0.8488 

0 . 7500 

0.6791 

0.6250 

0.5821 

0.5469 

1 

1.53^9 

1.0095 

.8065 

.6913 

.6143 

.5582 

.5U7 

.4797 

.4510 

2 

1.2500 

.8398 

.6797 

.5875 

.5250 

.4788 

.4428 

.4136 

.3895 

3 

1.0568 

.7208 

.5894 

.5130 

.4608 

.4218 

.3912 

•  3663 

.3453 

4 

.9167 

.6323 

.5213 

.4565 

.4ii8 

•  3784 

.3519 

•  3301 

.3118 

5 

.8102 

.5639 

.4679 

.4118 

•3731 

.3439 

.3206 

•  3014 

.2851 

6 

.7266 

.5092 

.4249 

.3756 

.3415 

.3156 

.2950 

.2779 

.2633 

7 

.6589 

.4645 

.389^ 

.3456 

.3152 

.2920 

.2735 

.2581 

.2450 

8 

.6031 

.4272 

.3596 

.3202 

.2928 

.2720 

•  2552 

.2413 

.2294 

00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4Note  that  a^s  and  f^s  are  given  for  both  even  and  odd  values  of  k 
and  s.  We  are  still  considering  only  symmetric  oscillations,  but  can  extend 
the  analysis  to  cases  where  k  and/or  s  are  odd  by  insisting  that  both  the 
static  moment  and  damping  moment  be  odd  functions  of  the  dependent  variable. 
This  is  done  by  employing  absolute  values  where  needed.  For  instance,  a  qua¬ 
dratic  static  moment  would  be  written  as  M  =  MQy  +  M^y |y  |  • 
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^ks 


\k 

sX 

0 

1 

2 

3 

4 

5 

6 

7 

8 

0 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1 

.9593 

.9575 

.95^2 

.9505 

.9468 

.9434 

•9399 

.9364 

•9337 

2 

.9375 

.9346 

.9297 

.9243 

.9188 

•9135 

•9084 

.9038 

.8997 

3 

.9247 

.9211 

.9155 

.9091 

.9027 

.8964 

.8905 

.8850 

•  8797 

4 

.9167 

•9127 

.9067 

.8998 

.8929 

.8861 

•8797 

•8735 

•8677 

5 

.9115 

.9073 

.9011 

.8940 

.8868 

.8798 

.8729 

.8665 

.8603 

6 

.9082 

.9039 

.8975 

.8904 

.8831 

.8758 

.8688 

.8621 

.8558 

7 

.9060 

.9016 

.8953 

.8882 

.8808 

•8735 

.8664 

.8596 

.8530 

8 

.9047 

.9003 

.8940 

.8868 

.8795 

.8722 

.8650 

.8582 

.8516 

00 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

Comparison  for  a  Cubic  Moment 
For  a  cubic  moment,  equation  (82)  looks  as  follows 


H0e  mQ  +  O.625O1112  mG  +  0.6797m2 
H0  mQ  +  0.9375m2  m0  +  0.9297m2 


(84) 


The  accuracy  of  this  expression  can  be  demonstrated  in  several  ways.  First, 
we  will  compare  equation  (84)  with  the  exact  solution  (7l)  or  (74)  when  only 
linear  damping  is  present  (h2  =  h4  =  .  .  .  =  hn  =  0) .  In  this  way  we  can  also 
compare  an  approximate  formula  that  is  given  by  Murphy  in  reference  7 >  which, 
in  our  notation,  is 


Hoe  _  rnQ  +  (3A)m2 
mo  +  (9/6)m2 


For  the  sake  of  comparison,  we  write  the  corresponding  case  of  (84)  with 
fractions  as 


Hoe  _  mp  +  (^/8)mg 
H0  ”  m0  +  (15/16)102 


(86) 


Both  of  these  formulas  have  the  correct  limit  as  m2  00  (Hoe/H0  =  2/3) , 
but  when  mQ  =  1,  formula  (85)  does  not  quite  have  the  correct  behavior  as 
m2  -*  0,  yielding  Hq  /hq  =  1  -  (3/8)m2  +  0(m22)  instead  of  the  exact  value 
H0e/H0  =  1  -  (5/l6)m2  +  0(m22)  of  the  present  analysis.  For  the  stable-stable 
moment  the  approximate  formula  (86)  is  practically  indistinguishable  from  the 
exact  solution  when  plotted  on  the  same  graph,  and  (85)  is  also  a  good  approx¬ 
imation  for  this  case.  These  results  are  shown  in  figure  6(a).  The  main 
deviations  occur  in  the  stable-unstable  and  unstable- stable  moments  near  the 
singular  points,  which  we  illustrate  in  figures  6(b)  and  (c).  The  approximate 
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(c)  Unstable -stable  1-3  static  moment. 


Figure  6.-  Comparison  of  approximate  solutions 
with  exact  solution. 

formula  (86)  of  the  present  method  is  seen  to  agree  very  well  with  the  exact 
curve  except  in  the  immediate  vicinity  of  the  singular  points,  which  are 
m2  =  -1  for  the  stable- unstable  moment  and  m2  =  2  for  the  unstable -stable 
moment.  The  present  approximation  (86)  is  obviously  superior  to  the  approxi 
mate  formula  (85)  . 


The  next  check  that  we  can  make  on  the  accuracy  of  equation  (84)  is  to 
allow  h2  to  he  nonzero.  A  number  of  different  cases  were  investigated,  and 
one  example  for  each  type  of  linear-cubic  moment  is  shown  in  figure  7 . 

Figure  7(a)  shows  results  for  a  stable-stable  moment.  The  magnitude  of  the 
nonlinearities  being  considered  on  this  plot  can  be  visualized  by  observing 


(a)  Stable- stable  1-3  static  moment;  stable - 
unstable  or  unstable -stable  1-3  damping 
moment . 


(b)  Stable -unstable  1-3  static  moment;  stable- 
unstable  or  unstable- stable  1-3  damping 
moment . 


(c)  Unstable- stable  1-3  static  moment;  stable - 
unstable  or  unstable -stable  1-3  damping 
moment . 


Figure  7*-  Comparison  of  present  analysis  with 
numerical  integrations . 
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the  small-scale  plots  of  static  and  damping  moments  that  appear.  These  are 
accurate  representations  of  the  nonlinearities,  not  sketches,  and  indicate  a 
wide  departure  from  a  linear  system.  The  first  thing  to  note  in  this  figure  is 
that  equation  (84)  agrees  to  this  scale  exactly  with  equation  (7l)  or  (74), 
which  in  turn  agrees  exactly  with  numerical  integration  of  equation  (l) .  The 
second  thing  to  note  is  that  a  negative  h2/m2  (recall  that  the  parameter 
h2/m2  is  a  constant,  independent  of  amplitude,  for  a  given  oscillator)  can  be 
obtained  in  two  ways,  H0  negative  and  H2  positive  or  H0  positive  and  H2 
negative.  The  first  case  leads  to  a  limit  cycle,  the  second  at  a  high  enough 
amplitude  to  an  unstable  situation.  The  surprising  thing  is  that  H0e/H0  is 
the  same  for  both  of  these  cases,  a  fact  that  was  revealed  in  the  analysis  but 
is  not  intuitively  obvious . 

Also  indicated  in  figure  7(a)  is  the  asymptotic  behavior  of  H0e/H0  as 
given  by  equation  (75) •  It  can  be  shown  that  expanding  equation  (84)  for  large 
m2  (with  h  ,h6,  .  .  .  =  0)  yields  the  same  values  as  equation  (75) >  except 
that  the  coefficient  of  m0h2/m2  is  O.289  instead  of  the  more  exact  value 
O.296. 

Figures  7(h)  and  (c)  show  representative  cases  for  the  stable -unstable  and 
unstable- stable  static  moments.  Once  again  the  magnitude  of  the  nonlinearities 
being  covered  is  shown  by  the  small-scale  plots.  For  both  of  these  cases,  the 
accuracy  of  equation  (84)  is  not  as  good  as  it  was  for  the  stable-stable 
moment,  but  the  agreement  is  still  excellent  except  in  the  immediate  vicinity 
of  the  singular  points.  The  final  point  to  be  made  is  that  the  curves  for  the 
limit  amplitude,  (-h2)z,  presented  in  figure  5,  can  be  closely  duplicated  by 
equation  (84)  except,  again,  in  the  immediate  vicinity  of  the  singular  points. 


Comparison  for  a  1-9  Moment 

A  more  extreme  nonlinear  case  was  next  investigated  to  check  the  adequacy 
of  the  approximate  solution  (82) .  A  1-9  static  moment  and  a  1-3  damping 
moment  (M0,M8,H0,H2  -f  0)  were  considered.  For  this  case,  equation  (82) 
appears  as 

H0g  _  m0  +  0.30l6me  +  m0  +  0.3596m8  ^ 

Ho  mQ  +  0.9047me  m0  +  0.8940m8 

In  this  case,  the  closed-form  solution  (69)  involves  hyperelliptic  integrals 
and  cannot  be  evaluated  by  simple  means,  but  a  comparison  can  be  made  with 
results  of  numerical  integrations  of  equation  (l) .  These  comparisons  are  made 
in  figure  8  for  a  stable-stable  and  unstable -stable  static  moment.  (The 
stable -unstable  case  is  not  too  interesting  in  that  the  moment  is  essentially 
linear  almost  to  the  singular  point  and  then  suddenly  goes  unstable.)  Instead 
of  the  parameter  h2/m2,  which  was  pertinent  for  the  cubic  moment,  the  param¬ 
eter  that  remains  constant  for  a  given  oscillator  is  now  h23|h2| /mB*  The 
comparisons  show  that  the  approximate  solution  is  not  so  good  as  it  was  for  a 
linear-cubic  static  moment.  However,  for  the  stable-stable  case,  agreement  is 
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still  excellent.,  For  the  unstable -stable  case ,  agreement  is  satisfactory 
except  close  to  the  singular  point  (i.e.,  m8  =  5)  • 


(a)  Stable- stable  1-9  static  moment;  1-3  (b)  Unstable -stable  1-9  static  moment;  1-3 

damping  moment.  damping  moment. 

Figure  8.-  Comparison  of  present  analysis  with 
numerical  integrations. 


Comparison  for  a  1-3-5  Moment 

The  final  comparison  that  was  made  was  for  a  l-3“5  static  moment  in  con¬ 
junction  with  several  different  damping  moments,  the  most  nonlinear  of  which 
was  also  1-3-5*  For  this  case,  equation  (82)  appears  as 

H0e  m0  +  O.625O1112  +  0. 4583m4  mQ  +  0.6797m2  +  0.5213m4 

Ho  mG  +  0.9375^2  +  0.9l67m4  mQ  +  0.9297m2  +  0.9067m4 


0.7500mo  +  0 . 525O1112  +  0.4ll8m4 

- x - 3  h 

mQ  +  0.9l88m2  +  0.8929m4 


(88) 


Too  many  possibilities  exist  to  cover  this  case  in  general.  Hence,  a  particu¬ 
lar  I-3-5  static  moment  was  selected  as  being  very  nonlinear,  a  stable- stab  le- 
unstable  moment.  This  moment  is  shown  in  figure  9(a).  Three  different 
damping  moments  were  then  chosen,  consisting  of  the  following:  (l)  linear 
term  only,  (2)  linear-cubic,  and  (3)  linear-cubic-quintic .  The  three  damping 
moments  are  shown  in  figure  9 ("b )  -  Note  that  they  were  selected  so  that  the 
average  value  of  the  damping  moment  between  0°  and  20°  was  the  same  in  each 
case . 
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(a)  l-3“5  static  moment. 


(b)  Damping  moments. 


Figure  9*"  Moments  used  for  more  complex 
comparisons. 


yQ»deg 


Figure  10.-  Comparison  of  present  analysis  with 
numerical  integrations;  l-3"5  static  moment; 
various  damping  moments. 


Equation  (88)  is  compared  with 
numerical 'integrations  of  equation  (l) 
for  the  three  cases  in  figure  10.  To 
the  scale  plotted  and  for  the  range 
covered,  no  differences  existed. 

All  of  these  comparisons  that 
have  heen  presented  indicate  that 
equation  (82)  is  surprisingly  accurate 
over  a  wide  range  of  nonlinearities  in 
hoth  the  static  and  damping  moment. 
With  this  established,  we  will  next 
concern  ourselves  with  determination 
of  the  coefficients  Mg  and  Hg  from 
data . 


DETERMINATION  OF  NONLINEAR  PARAMETERS  FROM  DATA 


Formula  (82)  with  (83)  can  he  useful  in  determining  the  nonlinear  damping 
parameters  of  an  oscillator  from  a  set  of  observed  oscillations.  Let  us  assume 
that  a  given  oscillator  is  governed  by  the  differential  equation  (l) .  We  wish 
to  infer  from  observed  oscillations  the  appropriate  values  of  HQ,H2,H4,  .  .  . 
and  MQ,M2,M4,  .... 

We  assume  that  the  damping  is  small  so  that  h  =  H Q/v  «  1  (or  more 
appropriately ,  since  v  =  |^yan|1/2  is  actually  arbitrary,  we  could  also  say 
Ho/03  «  1,  where  w  is  the  frequency  and  is  proportional  to  v)  .  This  assump¬ 
tion  is  normally  satisfied  in  ballistic-range  testing-  In  addition,  we  assume 
that  the  frequency  and  maximum  amplitudes  for  each  half  cycle  of  the  motion 
can  be  accurately  determined  from  the  data-  The  effective  linear  damping  for 
each  half  cycle  is  then  computed  by 


Ho( 


=  4  w  In 


n 


n+i 


(89) 


Using  an  approximation  for  the  frequency  quadrature  (eq.  (67b)  )  that  is 
developed  in  reference  13 ,  we  can  write  an  expression  for  the  frequency  as  a 
function  of  the  parameters  MQ,M1,M2,  ...  as 


where 


w2  =  MQ  + 


AnMnyan 


n>o 


(90a) 


n 


h 

n  +  2 


(90b) 


This  formula  is  a  good  approximation  for  large  and  small  nonlinearities,  except 
in  the  vicinity  of  the  singular  points  of  static  instability.  The  understand¬ 
ing  and  use^of  equation  (90a)  makes  it  possible  to  extract  the  appropriate 
values  of  Mn  by  fitting  equation  (90a)  to  a  set  of  frequency  versus  amplitude 
data  that  has  been  measured  for  a  given  oscillator.  For  instance,  if  an 
oscillator  is  governed  by  a  linear-cubic  static  moment,  then  equation  (90) 
yields 

a)2  =  MQ  +  ^  M2ya2 


The  frequency  squared  then  gives  a  straight  line  when  plotted  against  ya2 . 

The  slope  of  the  curve  yields  the^value  of  M2  and  the  w2  intercept  yields 
Mq.  When  h  «  1,  we  note  that  Mn  =  Mn[l  +  0(h2)].  Hence,  for  small  damping, 
the  values  of  the  actual  static  parameters  (Mn)  are  determined.  The  cubic 
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moment  is  the  simplest  case;  for  more  general  cases  the  static  moment  polyno¬ 
mial  (eq.  (90a))  is  terminated  so  that  an  adequate  fit  of  the  frequency  versus 
maximum  amplitude  data  is  obtained  (ref.  13)  • 

A  similar  but  slightly  more  complicated  situation  exists  for  determining 
the  damping  parameters  Hn-  In  this  case ,  the  static  parameters  play  an  impor¬ 
tant  role ,  whereas  for  the  determination  of  the  static-moment  parameters  Mn, 
the  effect  of  the  damping  parameters  can  be  neglected.  The  formula  for  H0e 
as  a  function  of  H0^H1,H 2.9  *  *  •  can  be  written  as 


Hoe  =  BoHo  +  Y,  B“Hnyan 
n>o 

where  equation  (69)  is  used  to  determine  Bn  as 


or,  from  equations  (82)  and  (83),  approximately  as 


°°ns^sya 


1  s>o 


Bn  2  +  n  yi 
s>o 


B  M  v  ' 
HnsA1s*ya 


(91a) 


(911) 


(91c) 


Expression  (91a)  is  analogous  in  form  to  (90a).  In  this  case,  however,  the 
coefficients  Bn  are  functions  of  MQ,M1,M2,  .  .  .  and  ya,  whereas  the  coef¬ 
ficients  An  are  constants.  In  the  determination  of  the  damping,  however,  the 
static  parameters  can  be  regarded  as  already  determined  by  the  use  of  (90)  . 
Hence,  the  coefficients  Bn  are  functions  of  ya  insofar  as  determination  of 
the  damping  is  concerned. 

It  is  now  useful  to  divide  both  sides  of  equation  (91a)  by  BQ,  obtaining 


n>o 


5s 

Bo 


Hny 


n 


n^a 


(92) 
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Because  .  .  .  each  vary  with  ya  in  a  similar  manner,  the  ratios 

Bn/B0  vary  slowly  with  ya.  Equation  (92)  can  now  be  used  to  study  a  set  of 
damping  data  in  a  manner  analogous  to  studying  frequency  data  with 
equation  (90) . 


As  a  simple  example,  assume  that  a  set  of  data  for  a  given  oscillator  has 
been  obtained  and  that  the  static  moment,  completely  arbitrary,  has  been  deter¬ 
mined  with  the  use  of  formula  (90) .  Let  it  be  assumed  that  the  damping  moment 

linear  term  and  a  single  arbitrary  nonlinear  term  so 


can  be  described  by  a 
that  (92)  appears  as 


^oe 

B 


=  H, 


o 


-  jjs  IW1 


(93) 


where  n  is  not  known.  Choosing  a  particular  value  of  n  determines  Bn/B0. 
Then  one  can  plot  H0  /B0  versus  (Bn/B0)yan.  When  the  right  value  of  n  is 
chosen,  the  data  will6 fall  on  a  straight  line.  The  slope  of  the  line  yields 
Hn,  and  the  intercept  of  the  line  with  the  Hoe/B0  axis  yields  HQ.  This  will 
be  true  regardless  of  the  particular  form  of  the  static  moment. 


For  an  arbitrary  damping  moment  described  by  a  given  set  of  data  over  an 
amplitude  range,  there  may  be  several  or  many  combinations  of  the  parameters 

.  .  .  that  will  fit  the  data.  This  situation  also  occurs  in  deter¬ 
mining  the  static-moment  parameters  (ref.  13)*  One  must  often  settle  for  a 
member  of  a  class  of  moments  that  gives  a  good  fit.  Thus,  a  certain  amount  of 


experience  may  be  useful  in  analyzing 


1 . . I _ L 


048  12  16 

ya 

Figure  11.-  Form  of  data  to  reveal  type  of 
nonlinearity  in  damping. 


xe  data.  For  a  set  of  data  including 
small  angles  of  ya  as  well  as  large, 
one  should  plot  the  damping  data  as 
Hoe/B0  versus  (B2/B0)ya2  (°r  as  a  first 
try  against  ya2  since  B2/bq  varies 
slowly  with  ya) .  If  the  data  fall  on 
a  straight  line,  the  damping  moment  is 
a  cubic  as  described  by  (93)  with 
n  =  2.  If  the  data  deviate  from  a 
straight  line,  which  might  be  expected 
to  occur  at  larger  values  of  ya,  one 
can  estimate  the  magnitude  and  sign  of 
the  next  higher  order  damping  term 
needed  to  fit  the  data.  In  such  a  man¬ 
ner  the  appropriate  form  of  damping 
polynomial  tends  to  suggest  itself. 

The  more  data  that  are  available, 
the  better  the  damping  moment  can  be 
defined.  Finally  the  data  can  be 
fitted  by  some  curve -fitting  technique 
for  various  combinations  of  damping 
parameters.  The  combinations  that  fit 
the  data  best  are  the  appropriate  ones 
to  represent  the  damping  moment. 

In  figure  11  the  results  of  fig¬ 
ure  10  are  plotted  in  the  manner  sug¬ 
gested.  Results  for  the  linear-cubic 
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damping  moment  fall  on  a  straight  line  as  expected.  The  deviation  of  the 
I-3-5  damping  data  from  a  straight  line  indicates  that  the  damping  moment  is 
more  complicated  than  linear  cubic;  the  fact  that  the  results  curve  downward 
from  a  straight  line  indicates  that  the  next  higher  order  term  required  is  a 
destabilizing  term,  which  is,  in  fact,  the  case. 


CONCLUDING  REMARKS 


The  preceding  investigation  has  considered  the  effects  of  nonlinearities 
on  damped  oscillations.  It  was  found  that  this  investigation  could  be  con¬ 
ducted  analytically  in  two  ways.  With  large  damping  nonlinearities  and  small 
“effective"  static-moment  nonlinearities,  an  approximate  solution  could  be 
established  by  means  of  an  integral  equation.  On  the  other  hand,  when  the 
damping  is  small  but  the  static  nonlinearities  arbitrarily  large  (a  common  sit¬ 
uation  in  aerodynamics),  an  approximate  solution  could  be  obtained  by  means  of 
an  equivalent  first-order  differential  equation. 

When  the  damping  is  small  and  the  static  nonlinearities  large,  it  is  con¬ 
venient  to  study  the  characteristics  of  the  nonlinear  oscillations  by  means  of 
a  parameter  called  the  "effective  linear  damping."  This  parameter  is  readily 
obtainable  from  experiment  and  offers  a  means  of  deducing  the  nonlinear  damp¬ 
ing  characteristics  of  an  oscillator  from  a  set  of  observed  oscillations. 
Comparisons  with  some  exact  solutions  and  with  numerical  solutions  showed  that 
the  approximate  formula  derived  herein  for  the  effective  linear  damping  is 
extremely  accurate  over  a  wide  range  of  nonlinearities.  The  limit  cycle  of  a 
nonlinear  oscillation  occurs  when  the  "effective  linear  damping"  is  zero,  and 
the  limiting  amplitude  is  easy  to  obtain  for  various  combinations  of  the  non¬ 
linear  damping  and  static  moments.  The  concept  of  an  "effective  linear 
damping"  is  thus  very  useful  in  analyzing  nonlinear  damping. 

The  integral  equation  (36)  and  the  first-order  equation  (39)  are  believed 
to  be  novel  formulations  for  nonlinear  oscillation  problems.  It  is  possible 
that  further  investigation  of  these  equations  would  be  fruitful. 


Ames  Research  Center 
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APPENDIX  A 


MODIFIED  EQUATIONS  OF  KRYLOFF  AND  BOGOLIUBOFF 


Here  we  wish  to  modify  the  approximation  methods  of  Kryloff  and 
Bogoliuboff  (ref.  4)  for  nonlinear  oscillations  so  that  the  pertinent  equations 
fit  conveniently  into  the  framework  of  the  present  analysis.  We  shall  forego 
the  Fourier  series  expansions  of  the  usual  Kryloff -Bogoliuboff  analysis  and 
find  forms  of  the  amplitude  and  frequency  equations  that  are  amenable  to  appro¬ 
priate  transformations  that  in  turn  yield  a  basic  exact  solution  as  a  special 
case . 


Following  Kryloff  and  Bogoliuboff,  consider  a  nonlinear  equation  of  the 
following  form : 


2 

+  V^X 


(Al) 


where  f  is  an  arbitrary  function  of  x  and  dx/dt.  If  f  =  0,  then  the  solu¬ 
tion  is  given  by 

x  =  a  cos(vt  +  cp)  (A2) 


where  a  and  cp  are  arbitrary  constants.  For  the  general  equation  (Al),  we 
assume  equation  (A2)  to  be  a  solution  with  a  =  a(t)  and  cp  =  cp(t)  suitably 
determined  functions  of  t . 

Let  us  choose  a(t)  and  cp(t)  by  following  the  method  of  variation  of 
parameters  for  linear  systems.  The  derivative  of  equation  (A2)  is 

—  =  -fr  cos(vt  +  cp)  -  a  (v  +  sin(vt  +  cp)  (A3) 

dt  dt  \  dt/ 

Since  we  have  two  functions  to  determine,  a(t)  and  cp(t) ,  we  have  two  conditions 
at  our  disposal.  One,  of  course,  will  be  the  original  differential  equa¬ 
tion  (Al) .  The  other  we  choose  as 


—  cos(vt  +  cp)  -  a  sin(vt  +  cp)  =  0 
dt  dt 


so  that  equation  (A3)  now  reads  simply 

=  -av  sin(vt  +  cp) 
dt 


(A4) 


(A5) 
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just  as  for  the  case  f  =  0.  Differentiating  equation  (A5)  and  substituting 
into  the  original  equation  (Al) ,  we  get 

sin( vt  +  9)  +  a  I?  cos(vt  +  9)  =  -  f[a  cos(vt  +  cp) ,  -av  sin(vt  +  cp)  ] 
dt  dt  v 

(A6) 

Expressions  (AA)  and  (A6)  are  two  equations  for  da/dt  and  dcp/dt .  Solving 
for  da/dt  and  d9/dt,  we  obtain 

|f  =  i  f [a  cos ( vt  +  cp)>  -av  sin(vt  +  cp)]sin(vt  +  cp)  (A7a) 

=  -1-  f[a  cos( vt  +  9),  -av  sin(vt  +  9)]cos(vt  +  9)  (A7b) 

dt  av 

These  are  the  basic  equations  of  Kryloff  and  Bogoliuboff. 

Equations  (A7a,  b)  are  coupled  and  thus  must  be  solved  simultaneously. 
Equation  (A7a)  can  be  uncoupled,  however,  if  we  introduce,  in  place  of  t,  a 
new  independent  variable  t  defined  as  follows: 

T  =  vt  +  9  (A8a) 

It  follows  that 

±r  =  _ v -  (A8b) 

dt  dcp 

dT 


Equations  (A7a,  b)  may  now  be  written  in  terms  of  r  as  follows 


da 

dT 


i  f (a  cos  t ;  -av  sin  r)sin  t 
v2 


1  +  — f (a  cos  t y  -av  sin  t)cos  t 
av2 


(A9a) 


— f(a  cos  t,  -av  sin  t)cos  t 
dcp  _  av2  _ 

dT  n  1 


1  +  — -  f (a  cos  t y  -av  sin  t)cos  t 


av 


(A9b) 


=  v  1  +  -i-  f (a  cos  t  .  -av  sin  t)cos  t 
dt  av2 


(A9c) 


We  shall  term  equations  (A9a,  b,  c)  the  modified  Kryloff -Bogoliuboff  equa¬ 
tions.  Equations  (A9a)  and  (A9c)  are  the  pertinent  ones  since  the  solution  to 
equation  (Al)  is  given  by  y  =  a(r)cos  t-  Equation  (A9a)  is  a  first-order 
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nonlinear  equation  for  a  as  a  function  of  t,  and  it  is  uncoupled  from  the 
other  two.  When  it  has  been  solved,  the  other  equations  can  be  solved  at  least 
by  quadratures. 
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